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Abstract 

This  paper  presents  a  lumped  element  model  of  a  natural  gas  steam  reformer  being  part  of  a  5  kW(el.)  residential  fuel  cell  power  plant.  It 
describes  the  thermal  behaviour  over  the  whole  operating  range  including  plant  startup.  For  the  reforming  reactions  chemical  equilibrium  is 
assumed.  The  derivation  of  the  model  starting  with  basic  thermal  and  chemical  laws  is  shown.  Therefore  model  parameters  have  a  physical 
meaning  and  were  calculated  with  the  assumption  of  lumped  system  properties.  To  achieve  a  better  agreement  with  experiments  the  parameters 
were  furthermore  identified  from  measured  data.  Experimental  data  of  a  demonstration  plant  were  used  to  verify  the  model. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Fuel  cell  power  plants  are  a  promising  technology:  Using 
hydrogen  as  fuel  they  provide  a  local  emission  free  energy 
source  with  high  electrical  efficiency.  To  increase  overall  ef¬ 
ficiency  also  the  waste  heat  should  be  used  in  a  combined 
heat  and  power  plant.  Therefore  currently  small  power  plants 
for  residential  energy  supply  are  developed  which  should  re¬ 
place  the  conventional  heating  system.  But  due  to  the  lack  of 
a  hydrogen  infrastructure  the  already  available  natural  gas  is 
the  fuel  for  the  near  future.  If  using  natural  gas,  hydrogen  is 
produced  by  a  reforming  process. 

Recent  developments  show  the  feasibility  of  such  small 
units  and  some  demonstration  plants  already  exist.  Up  to  now 
only  very  limited  data  of  these  plants,  including  the  reforming 
process,  are  published. 

In  this  paper  a  dynamic  model  of  a  steam  reformer  for  a 
residential  fuel  cell  system  is  presented.  The  model  was  con¬ 
structed  using  basic  physical  and  chemical  laws  and  therefore 
provides  a  good  insight  into  the  process  and  covers  the  com¬ 
plete  operating  range  including  startup  and  shutdown.  This 
considerations  may  be  used  also  for  different  reformer  types. 
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Our  applications  are  control  system  design  and  the  opti¬ 
mization  of  the  plant  startup  procedure.  The  first  application 
can  be  divided  into  controller  design  and  observer  design:  For 
controller  design  the  model  is  only  applied  offline  during  the 
design  process.  In  difference,  in  an  observer  the  model  is  ap¬ 
plied  online,  i.e.  it  runs  in  parallel  to  the  real  plant  and  is  used 
to  estimate  plant  data.  This  may  be  the  outlet  gas  composition 
which  can  not  be  measured  directly  with  acceptable  expense 
under  real  operating  conditions.  The  second  application,  op¬ 
timization  of  plant  startup,  is  targeted  on  the  calculation  of 
control  inputs  that  minimize  some  objective  function.  There¬ 
fore  a  model  of  the  whole  plant,  with  the  reformer  as  one 
component,  is  required. 

In  the  first  part  of  this  paper  the  demonstration  plant  is 
described.  In  the  second  part  the  model  structure,  assump¬ 
tions  and  simplifications  are  discussed.  The  third  part  covers 
the  theoretical  calculation  and  identification  of  the  model  pa¬ 
rameters  and  finally  the  model  is  compared  with  real  plant 
data. 


2.  Demonstration  plant 

The  described  reformer  is  shown  in  Fig.  1.  It  is  part  of  a 
demonstration  plant  of  a  residential  fuel  cell  power  station 
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Nomenclature 

C 

heat  capacity  (c  •  m) 

c 

specific  heat  capacity 

G 

Gibbs  enthalpy 

g 

specific  Gibbs  enthalpy 

Ah 

heat  of  reaction 

k 

factor  for  heat  conduction/radiation 

m 

mass 

n 

amount  of  substance 

An 

extent  of  reaction 

h 

molar  flow  rate 

P 

pressure 

P0 

standard  pressure 

P 

parameter  vector 

Q 

heat  flow 

R 

gas  constant 

r 

heat  of  evaporation 

S/C 

steam  to  carbon  ratio 

s 

valve  position 

T 

temperature 

u 

parameter  vector 

y 

concentration 

X 

excess  air  value 

Indices 

A 

ambient 

B 

burner 

c 

convection 

E 

evaporator 

e 

evaporation 

F 

flame 

G 

ground  plate 

i 

inlet 

0 

outlet 

R 

reactor 

r 

reaction 

V 

burner  valve 

w 

wall 

located  at  the  Fachhochschule  Ulm.  The  plant  was  devel¬ 
oped  by  the  Centre  for  Solar  Energy  and  Hydrogen  Research 
(ZSW),  Ulm,  and  the  gas  processing  subsystem  (described 
in  [1])  by  the  Fraunhofer  Institute  for  Solar  Energy  Systems 
(ISE),  Freiburg. 

The  plant  components  are  shown  in  Fig.  2.  Natural  gas 
and  water  are  fed  to  the  reformer  which  is  heated  by  a  natural 
gas  burner.  The  product  gas  is  a  mixture  containing  mainly 
hydrogen,  carbon  dioxide,  some  methane  and  also  carbon 
monoxide  in  the  range  of  a  few  percent.  As  carbon  monoxide 
acts  as  a  catalyst  poison  for  the  fuel  cell,  it  is  removed  in 
a  gas  purification  system  consisting  of  two  shift  stages  and 
a  PROX-unit.  Before  entering  the  fuel  cell  stack  the  gas  is 
humidified.  The  stack  is  a  40  cell  PEM  type  with  5  kW(el.) 


Fig.  1.  The  steam  reformer.  Approximate  dimensions  are  60  cm  height  and 
28  cm  diameter. 

nominal  power.  The  hydrogen  which  is  not  consumed  in  the 
cell  is  fed  to  a  catalytic  burner. 

A  detailed  cross-section  of  the  reformer  is  depicted  in  Fig. 
3.  Natural  gas  (methane)  and  liquid  water  enter  the  evapora¬ 
tor  section  at  the  top  and  flow  to  the  center  inside  a  helical 
tube.  The  gas  and  water  are  heated  by  the  hot  flue  gas  from 
the  gas  burner.  The  methane/steam  mixture  then  enters  the 
reactor  consisting  of  an  inner  and  outer  tube  filled  with  cata¬ 
lyst  material.  At  the  catalyst  the  following  two  reactions  take 
place  at  temperatures  above  500  °C  [2]: 

CH4  +  H20  ^  CO  +  3H2  Ah0  =  206  kJ/mol  (1) 

CO  +  H20  ^  C02  +  H2  A/zi  =  -41  kJ/mol  (2) 

The  outlet  composition  depends  on  the  operating  conditions 
and  is  typically  about  10%  C02,  <10%  CO,  a  few  percent 
CH4  and  the  remainder  H2  (dry  gas).  The  overall  reaction 
is  endothermic.  The  reaction  heat  is  provided  by  a  natural 
gas  burner  located  around  the  reactor  tubes.  Natural  gas  and 
air  are  mixed  outside  the  burner  and  flow  through  perforated 
ceramic  plates.  The  combustion  takes  place  at  the  inner  plate 
surface: 

CH4  +  202  — >  C02  +  2H20  A hc  =  -804  kJ/mol 

(3) 

The  flue  gas  flows  through  the  gap  between  the  two  reactor 
tubes  and  then  through  the  evaporator.  Therefore  the  reactor 
is  heated  by  radiation  from  the  burner  plates  and  convection 
of  the  flue  gas. 

The  inlet  (CH4)  gas  flow  rate  can  be  varied  between  10 
and  25  SLPM  (standard  liter  per  minute).  The  water  flow 
rate  is  given  through  the  molar  ratio  of  water  to  methane 
(steam  to  carbon,  S/C)  which  is  approximately  3.5.  Tem¬ 
peratures  are  measured  at  the  evaporator  outlet,  the  reac¬ 
tor  outlet  and  at  four  different  locations  of  the  burner  using 
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Fig.  2.  Components  of  the  demonstration  plant. 


K-type  thermocouples.  The  outlet  gas  composition  is  ana¬ 
lyzed  with  an  Emerson  NGA2000  process  gas  analyzer  (con¬ 
taining  NDIR-(infrared-)sensors  for  CO2,  CO,  CH4  and  a 
heat  conduction  sensor  for  H2). 

The  plant  is  controlled  by  an  industrial  programmable 
logic  controller  (PLC).  For  the  reformer  system  it  calculates 
the  natural  gas  and  water  flow  rates  and  controls  the  burner 
temperature.  The  reforming  gas  flow  rate  is  set  with  a  speed 
controlled  compressor  and  the  water  flow  rate  with  a  dosing 
pump. 

The  experimental  data  comprise  two  different  operation 
points,  i.e.  40  and  70  SLPM  reformate  flow,  each  with  step¬ 
wise  variations  in  gas  and  water  flow  rates  and  burner  tem¬ 
perature.  The  operating  pressure  is  3.5  bar. 

3.  Reformer  model 

Inputs  and  outputs  of  our  model  are  listed  in  Table  1. 
Before  continuing  with  its  description  some  common  ap- 


inlet  (methane,  water) 


Fig.  3.  Cross-section  of  the  reformer. 


proaches  shall  be  shortly  reviewed.  Models  of  dynamic  sys¬ 
tems  can  be  divided  into  theoretical  and  experimental  ones 
[3].  Theoretical  models  rely  on  physical  equations  and  there¬ 
fore  they  describe  the  inner  structure  of  a  system.  Usually 
they  are  expressed  as  differential  equations. 

In  contrast  experimental  models  are  based  on  measured 
data  (e.g.  the  step  response)  and  therefore  they  describe  the 
relation  between  inputs  and  outputs,  but  not  the  inner  struc¬ 
ture  of  a  system.  Based  on  the  experimental  data  the  parame¬ 
ters  of  some  standard  structure  (e.g.  a  transfer  function)  may 
be  calculated.  This  is  called  fitting  or  identification  and  the 
result  is  a  parametric  model  [4]. 

In  practice  often  some  intermediate  approaches  are  used, 
for  example,  the  structure  of  a  parametric  model  results  from 
theoretical  modeling  but  its  parameters  are  identified  based 
on  experimental  data. 

Examples  are  Meziou  and  Alatiqi  [5],  which  describe  the 
development  of  a  control  system  for  an  industrial  steam  re¬ 
former.  First  they  give  a  theoretical  model  of  this  reformer  and 
identify  its  parameters  based  on  experimental  data.  Then,  for 
control  system  design,  they  derive  a  linear  ARX  model  from 
the  initial  model  (ARX  is  a  standard  structure  consisting  of  a 
discrete  time,  rational  transfer  function  and  delay  time).  Such 
a  model  is  well  suited  for  control  system  design  if  operating 
conditions  are  close  to  the  one  used  during  identification,  but 
we  also  need  to  describe  the  reformer  startup. 

An  example  of  a  theoretical  steam  reformer  model  can 
be  found  in  [6].  This  paper  describes  a  tubular  fixed  bed 
reactor  and  considers  reaction  kinetics,  diffusion  and  heat 
conduction.  This  results  in  a  two-dimensional  partial  differ¬ 
ential  equation  (PDE).  For  simulation  a  finite-difference  dis¬ 
cretization  is  applied.  This  model  only  includes  the  reactor 
itself  but  no  further  components  like  evaporator  and  burner. 
A  one-dimensional  model  of  an  auto-thermal  methanol  re¬ 
former  for  automotive  use  is  described  by  Maume  [7].  We 
do  not  use  a  PDE  structure  because  of  the  much  more  com¬ 
plex  geometry  of  our  system  and  because  simple  ordinary 
differential  equation  (ODE)  models  are  preferred  for  control 
system  design. 

General  information  on  the  chemistry  of  steam  reforming 
can  be  found  in  [8]  and  targeted  to  fuel  cell  systems  in  [2]. 
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Table  1 


Model  inputs  and  outputs 


Inputs 

— 

Burner  state  (fan  on/off,  gas  valve  open/closed) 

Sy 

Burner  valve  position 

nn2o,  Th2o 

Water  flow  rate  and  temperature 

^ch4  ,  Tch4 

Natural  gas  flow  rate  and  temperature 

Outputs 

7w 

Wall  temperature 

Tb 

Burner  temperature 

Te 

Evaporator  temperature 

Tr 

Reactor  temperature 

nR 

Outlet  gas  flow  (CH4,  H2O,  H2,  CO2,  CO) 

The  physical  laws  are  generally  represented  as  partial  dif¬ 
ferential  equations.  But  to  ensure  later  applicability  for  con¬ 
trol  system  design  and  also  as  observer  in  an  operational 
system  a  lumped  model  structure  in  the  form  of  ODEs  is 
used. 

The  lumped  elements  are  wall,  ground  plate,  burner,  re¬ 
actor  and  evaporator.  Each  of  them  is  assigned  one  uniform 
temperature: 

•  wall  (temperature  7w):  This  is  the  tube  around  the  burner. 
The  corresponding  temperature  sensor  is  mounted  be¬ 
tween  the  wall  and  the  burner  and  indeed  measures  the 
temperature  of  the  burner  gas,  which  is  assumed  to  adopt 
the  wall  temperature. 

•  ground  plate  (7g):  The  steel  plate  at  the  bottom  of  the 
reformer  assembly.  Its  temperature  is  not  measured. 

•  burner  (7b):  The  burner  consists  of  ceramic  plates  located 
around  the  reactor.  The  temperature  sensors  are  near  the 
plate  surface  where  combustion  takes  place.  It  is  assumed 
that  the  sensors  correspond  to  the  surface  temperature. 

•  evaporator  (7p):  It  is  the  top  part  of  the  reformer.  Here  the 
outlet  temperature  (equal  to  the  reactor  inlet  temperature) 
is  of  interest. 


•  reactor  (7r):  It  is  located  in  the  center  of  the  burner. 
There  is  no  temperature  sensor  inside  the  reactor,  but  as 
the  gas  composition  depends  on  the  reactor  temperature 
(discussed  below),  the  reactor  temperature  is  defined  by 
the  corresponding  outlet  gas  composition.  A  temperature 
sensor  located  at  the  reformer  outlet  only  approximately 
agrees  with  the  actual  reactor  temperature,  especially  not 
during  dynamic  operation  (e.g.  plant  startup). 

In  the  following  the  heat  exchange  mechanisms  between 

the  above  elements  are  discussed  (see  also  Fig.  4): 

•  Heat  conduction  is  modeled  as  proportional  to  the  temper¬ 
ature  difference  with  a  proportionality  factor  k.  k  results 
from  geometry  and  heat  conductivity  and  is  constant  (i.e. 
not  temperature  dependent).  The  heat  transfer  from  the 
evaporator  to  the  ambience  consists  of  conduction  in  the 
thermal  insulation  and  free  convection.  These  two  effects 
are  not  modeled  separately. 

•  Radiation  takes  place  between  burner  and  reactor.  The  heat 
flow  is  proportional  to  the  difference  of  the  4.  power  of 
the  temperatures  with  the  proportionality  factor  &br  de¬ 
scribing  geometry  and  surface  emissivity  (which  is  as¬ 
sumed  to  be  temperature  independent).  Radiation  from 
burner  to  wall  is  negligible,  because  the  outside  sur¬ 
face  temperature  of  the  burner  plates  is  much  lower 
than  their  inside  temperature.  The  temperature  gradi¬ 
ent  over  the  burner  cross-section  [9]  is  not  taken  into 
account. 

•  Convective  heat  transfer  is  considered  for  both  burner  and 
reformate  gas  streams.  But  its  calculation  is  only  feasi¬ 
ble  for  simple  geometries  [10].  Therefore  and  to  keep  the 
model  simple  it  is  assumed  that  gas  flowing  through  or 
along  an  element  adopts  the  element  temperature.  Then 
the  heat  flow  can  be  calculated  from  an  energy  balance. 
There  is  only  one  exception:  It  is  not  likely  that  the  flue 
gas  cools  down  to  the  ground  plate  temperature.  Therefore 
here  a  constant  kpG  is  introduced,  which  corresponds  to 


Fig.  4.  The  model  structure:  The  boxes  correspond  to  the  lumped  elements,  single  lines  depict  heat  transfer  (solid:  conduction,  dashed:  radiation),  double  lines 
depict  the  burner  gas  flow  and  triple  lines  the  reformate  gas  flow. 
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the  fraction  of  heat  transferred.  So  the  flue  gas  cools  down 
to  Tpc  =  Tp  —  &fg(Tf  —  7g). 

•  The  combustion  takes  place  at  the  surface  of  the  burner 
which  is  heated  by  the  hot  flue  gas.  This  heat  flow  is  as¬ 
sumed  to  be  proportional  to  the  difference  of  flame  and 
burner  temperature. 

•  Methane  and  water  (liquid)  entering  the  evaporator  are 
heated  up  to  the  evaporator  temperature,  i.e.  the  water 
evaporates. 

•  The  methane  steam  mixture  leaving  the  evaporator  with 
7e  enters  the  reactor  and  is  heated  up  to  the  reactor  tem¬ 
perature  Tr.  The  reaction  takes  place  at  this  temperature. 

There  are  two  approaches  for  the  calculation  of  the  outlet 
gas  composition.  The  first  is  based  on  reaction  kinetics  (used 
in  [6,7,1 1])  and  requires  the  solution  of  a  ODE  system,  while 
the  second  assumes  chemical  equilibrium  ([12])  and  leads 
to  a  nonlinear  algebraic  equation.  We  use  the  second  method 
and  achieve  reasonable  results  for  the  flow  rates  applied  in  our 
system.  Of  course  at  higher  flow  rates  the  influence  of  reaction 
kinetics  would  become  apparent,  but  this  is  not  within  the 
operation  range  of  our  system. 

The  equilibrium  composition  is  defined  by  minimal  free 
enthalpy.  Its  practical  calculation  is  explained  briefly.  Be¬ 
cause  no  analytical  solution  exists  and  the  numerical  solution 
is  computationally  expensive,  also  an  approximate  solution  is 
presented.  This  is  better  suited  for  online  application  as  part 
of  a  plant  control  system  with  limited  computation  resources 
and  real-time  requirements. 

Consider  a  gas  mixture  with  all  N  components  of  the  re¬ 
forming  reactions  (1,  2)  at  constant  pressure  p ,  temperature 
7r  and  initial  composition  ncn4,i  and  mh2o,/  (Natural  gas 
may  contain  considerable  amounts  of  other  gases  including 
hydrocarbons  other  than  methane).  For  our  gas  supply  this 
is  not  the  case.  Otherwise  these  further  components  must  be 
included  in  the  free  enthalpy  calculation  and  maybe  other  re¬ 
actions  than  (1)  and  (2)  can  take  place).  As  the  two  reactions 
are  independent,  all  possible  compositions  can  be  expressed 
using  the  extent  of  reaction  A  no  and  A  n\\ 


Hr  = 


/  «ch4,r  \ 

(  ncu4,i  \ 

(-1 \ 

«h2o,r 

nu20  ,i 

-1 

”h2.r 

— 

0 

+  Amo 

3 

«co2,r 

0 

0 

\  ^CO,R  / 

0  / 

l  1  / 

+An\ 


0  \ 

-1 

1 

1 

K-i/ 


(4) 


Here  n  denotes  the  amount  of  substance.  The  range  of  Amo 
and  An i  is  restricted  by  the  condition  m  >  o.  With  these 


definitions  the  free  enthalpy  G  can  be  written  as  ([13]): 


N 


G{A«0,  A«i}  =  (  gi{T }  +  RT In 


i=  1 


Pi 

P0 


Hi,  R 

Pi  =  Ps-N  „ 

2^=1  nU  R 

The  equilibrium  corresponds  with  the  minimum  of 
G{Amo,  Ami}  which  is  denoted  by  G*: 


G*{Amq,  Am^}  =  min  G  (5) 

A^o,A^l 

The  minimum  is  determined  numerically  using  an  opti¬ 
mization  method  based  on  a  random  search  algorithm  with 
initial  conditions  Amo  =  Ami  =  0-  The  transition  to  molar 
flows  is  straightforward  by  replacing  n  by  fi.  The  equilibrium 
calculation  was  verified  with  data  from  [2]  and  [8]. 

Also  a  polynomial  approximation  for  the  equilibrium  with 
a  structure  based  on  the  two-dimensional  Taylor  series  was 
derived.  Operation  at  constant  pressure  is  assumed  and  there¬ 
fore  the  pressure  dependence  is  ignored.  The  approximation 
{Amq,  A h\]  is  only  valid  within  a  certain  temperature  and 
S/C-range: 


A"o  =  P0  ■  “  •  «H20,i 
a"i  =  Pi  •  “  •  «h2o ,i 
with 


u  =  (1  a  a2  a 3  /3  /32  ft 3  aft  a2/3  ap2)T 

(7) 


a  = 


Tr 

lOOtf 


^H2Q  ,i 

ncH4,i 


-3.5  = 


5 

C 


-3.5 


(8) 

(9) 


Scaling  of  S/C  and  Tr  is  optional  but  improves  numeri¬ 
cal  accuracy.  The  parameters  and  p  were  fitted  with 
the  least  squares  method  for  p  =  3 bar,  S/C  =  [2;  5]  and 
Tr  =  [500  °C;  800  °C]: 


po  =  (195.0  88.22  -5.504  -9.538  -30.41  7.821 

-2.223  -27.16  -4.443  7.684)  •  1(T3  (10) 


£j  =  (134.5  14.02  -19.62  2.491 

-11.94  0.09909  0.3631  0.7817  2.711 

—2.1 10)  •  10-3  (11) 

Fig.  5  contains  the  exact  (Eq.  (5))  and  approximated 
equilibrium  concentrations  for  p  =  3 bar  and  S/C  =  3. 
Deviations  over  the  whole  Tr  and  S/C- range  are  widely  less 
than  1%,  the  maximal  deviation  is  2.2%  for  CH4. 

Fluid  dynamics,  i.e.  transport  delay  and  pressure  drop  are 
not  included  in  the  presented  model,  because  these  are  small 
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Fig.  5.  Equilibrium  concentrations  for  p  =  3  bar  and  S/C  =  3  (dry  gas),  the  solid  lines  are  calculated  based  on  the  minimum  of  the  free  enthalpy  and  the 
dashed  lines  are  a  polynomial  approximation. 


compared  to  the  following  plant  components  and  therefore 
can  be  neglected. 

The  description  above  results  in  the  following  differential 
equations  taking  into  account  the  heat  capacities  C,  the  heat 
exchange  coefficients  k ,  the  burner  gas  flow  rate  Ab,  the  flue 
gas  flow  rate  Ap  and  the  reformer  gas  flow  rate  Ar  : 

wall: 

Cw  —  Tw  =  £gw(7g  ?w)  ^wa(Tw  7a) 

at 

~cp,b  •  «b(7w  -  7a)  (12) 

ground  plate: 
d 

Cg-7g  =  £bg(7b  7q)  ^gwCTg  7\y) 

at 

+/:fg  •  cPt f  •  hp(Tp  —  Tq)  (13) 

burner: 

Cb  — 7b  =  ^fb(7f  —  7b)  —  cp,b  •  Ab(7b  —  7\y) 

-*br(7b  -  7r)  -  *BG(rB  -  7g)  (14) 

Table  2 

Burner  and  flue  gas  flow  rates  depending  on  the  burner  state 


Component 

Burner  off 

Fan  on 

Burner  on 

«b,o2 

0 

2Any 

2Xh\ 

f*B,N2 

0 

2  A.fly 

2  yy  Any 

^b,ch4 

0 

0 

fly 

hpo2 

0 

2Xh\ 

2(A  —  l)fly 

tip  n2 

0 

2  Yj  A-fly 

2  yy  Any 

tl  pco2 

0 

0 

fly 

tip  h2o 

0 

0 

2h\ 

evaporator: 

d 

Ce  —  Te  =  £re(7r  -  Tfl)  +  cPj f  •  «f(7r  —  Tfl) 

-r  •  Ah2o,/  -  <a?,h2o  •  ^h2o(7e  -  7h2o) 

—Cp, ch4  •  ^ch4(7e  -  7ch4) 

-*ea(7e  -  7a)  (15) 

reactor: 

Cr  — 7r  =  kBR(Tg  -  Tr)  —  £re(7r  —  Tfl) 

+cp,f  •  ^f(7fg  —  7r)  —  AAo  •  AAo 
—Ah\  •  AAi  -  cp?h2o  •  ^h2o(7r  -  Tfl) 
-<a?,ch4  •  ^ch4(7r  -  7e)  (16) 

Table  3 


Model  parameters 


Parameter 

Theoretical  Value 

Identified  Value 

Unit 

Cw 

9.2  x  103 

7.27  x  103 

J/K 

CB 

0.56  x  103 

0.22  x  103 

J/K 

Ce 

9.35  x  103 

5.42  x  103 

J/K 

Cg 

3.5  x  103 

2.44  x  103 

J/K 

Cr 

4.65  x  103 

3.61  x  103 

J/K 

^0 

— 

-2.46  x  10“4 

mol/s 

iq 

— 

4.10  x  10“3 

mol  V/s 

A 

— 

1.7 

— 

AbR 

(0.78... 2.2)  x  10“9 

1.32  x  10“9 

W/K4 

AbG 

0.83...  2.8 

4.50 

W/K 

Agw 

0.77...  2.3 

5.16 

W/K 

Aea 

0.11. ..2.7 

0.439 

W/K 

Are 

1.1. ..4.3 

16.3 

W/K 

Awa 

0.32...  6.4 

1.16 

W/K 

Afb 

10 

16.1 

W/K 

AfG 

0.5 

0.30 

— 
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with 


Tp  =  Tq 


A/zc  •  ^b,ch4 
Cp, CH4  •  +  ^FB 


(17) 


resulting  from  the  energy  balance  of  the  combustion. 


valve  position.  If  the  burner  is  in  operation,  the  right  column 
is  valid.  The  flue  gas  flow  is  calculated  via  the  combustion  re¬ 
action  (3).  Other  than  the  above  listed  three  operation  modes 
do  not  occur. 

The  model  is  implemented  in  Matlab/Simulink. 


The  burner  gas  flow  rate  hB  and  the  flue  gas  flow  rate  hF 
depend  on  the  burner  valve  position  sy  and  the  burner  state 
(fan  and  gas  shut-off  valve)  and  are  listed  in  Table  2.  Here  vq 
and  v\  describe  the  linearized  valve  characteristics: 

n\  =  v  o  +  v\ sy  (18) 

Obviously  if  the  burner  is  off  (left  column)  the  flow  rate  is 
zero.  If  only  the  fan  is  on  (this  is  the  case  at  burner  startup  and 
shutdown)  the  air  flow  depends  on  the  X  setting  and  the  control 


4.  Model  parameters 

Table  3  depicts  all  parameters  with  the  theoretically  cal¬ 
culated  and  identified  (measured)  values. 

At  first  the  parameters  were  theoretically  determined.  For 
some  of  them  a  minimal  and  a  maximal  value  is  given  to 
account  on  uncertainty.  For  example  consider  the  radiation 
from  the  burner  to  the  reactor.  The  geometry  of  the  burner 


a) 

valve  position  sy 

b) 

* 

inlet  gas  flow  rate  UcH4,i 

0 

inlet  water  flow  rate  Vo,* 

c) 

* 

wall  temperature  TV 

0 

burner  temperature  Tb 

d) 

* 

evaporator  temperature  If 

0 

reactor  temperature  If 

e) 

* 

H2  concentration 

0 

CO 2  concentration 

f) 

* 

CO  concentration 

0 

CH4  concentration 

Fig.  6.  Simulation  with  operating  point  at  40%  of  nominal  power  (solid  lines  are  experimental  data,  dashed  lines  are  simulated). 
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and  the  reactor  are  not  as  simple  that  a  precise  determination 
of  &br  is  feasible.  Moreover  the  effective  emissivity  could 
only  be  estimated  within  a  certain  range.  The  parameters  of 
the  burner  valve,  vo  and  v\  are  determined  directly  from  mea¬ 
sured  burner  gas  flow  rates  at  different  valve  positions.  The 
A  value  is  calculated  from  the  flue  gas  composition. 

As  simulations  with  these  theoretical  values  did  not  show 
a  perfect  agreement  with  measured  data,  parameters  were 
subsequently  identified  for  the  predefined  model  structure. 
The  theoretical  parameters  were  used  as  initial  values  for  the 
identification  algorithm  fmincon[14].  The  objective  was  to 
minimize  the  quadratic  difference  between  simulation  out¬ 
put  and  experimental  data  (consisting  of  7b,  7e,  7w,  ycH4 
and  yco  with  appropriate  scaling  factors).  We  used  an  ex¬ 
perimental  data  set  comprising  four  sequences  with  different 
operating  conditions. 

5.  Results 

To  demonstrate  the  model  reliability,  simulated  and  mea¬ 
sured  behaviour  were  compared.  Figs.  6  and  7  show  the  sim¬ 


ulation  results  of  two  different  experiments.  Both  have  a  sim¬ 
ilar  sequence:  In  diagrams  (a)  and  (b)  the  model  inputs  are 
shown.  These  are  the  burner  control  valve  position  (a)  and  the 
water  and  natural  gas  flow  rates  (b).  All  further  diagrams  show 
the  model  outputs  (dashed)  in  comparison  to  the  respective 
measured  values  (solid).  All  relevant  temperatures  are  in  (c) 
and  (d).  The  reactor  temperature  is  not  measured,  but  instead 
the  reactor  outlet  temperature  is  shown.  Note  that  the  outlet 
temperature  does  only  approximately  and  only  during  near 
steady- state  operation  match  the  actual  reactor  temperature. 
It  is  not  meaningful  during  startup.  The  last  two  diagrams 
(e)  and  (f)  contain  the  gas  composition  on  a  dry  basis.  As 
the  concentrations  depend  on  the  reactor  temperature  they 
indirectly  confirm  the  thermal  model  behaviour. 

The  first  experiment  (6)  contains  the  reformer  startup  (at 
about  70  min)  and  operation  at  40%  of  the  nominal  power. 
One  after  another  the  natural  gas  flow  rate  (at  120  min),  the 
water  flow  rate  (at  160  min)  and  the  burner  temperature  were 
increased  (at  230  min).  The  second  experiment  is  similar,  but 
at  70%  of  nominal  power  and  with  decreased  flow  rates. 

The  temperatures  generally  show  a  good  agreement.  Dif¬ 
ferences  in  the  burner  temperature  (c)  are  hardly  to  notice. 


0.04 
.to 

|  0.03 

CD 

as  0.02 
£ 

§L  0.01 
0 

0  100  200  300  400  500 

(b)  t/min. 


Fig.  7.  Simulation  with  operating  point  at  70%  of  nominal  power.  See  Fig.  6  for  the  legend. 
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Indeed  during  startup  the  simulated  value  is  somewhat  to 
high  while  for  steady  state  conditions  it  is  somewhat  to 
low.  At  startup  differences  can  be  seen  at  the  wall  tem¬ 
perature  (c):  in  reality  this  temperature  rises  with  an  obvi¬ 
ous  delay.  The  difference  is  caused  by  the  lumped  structure 
which  cannot  take  into  account  temperature  inhomogene¬ 
ity  during  startup.  As  the  wall  temperature  is  not  of  ma¬ 
jor  relevance  for  the  system  behaviour  this  deviation  can  be 
tolerated. 

Diagram  (d)  depicts  the  evaporator  temperature.  The  in¬ 
fluence  of  changes  of  the  water  flow  rate  is  smaller  in  the 
model  than  in  reality.  This  can  be  explained  by  the  re¬ 
duction  of  the  heat  exchanger  to  one  lumped  element.  But 
a  more  sophisticated  model  seems  not  necessary  because 
this  temperature  is  not  an  important  value  for  the  plant 
control. 

The  gas  concentrations  are  also  matched  well  (diagrams  e 
and  f),  but  obvious  differences  exist  for  methane.This  is  due  to 
the  high  temperature  sensitivity  of  the  chemical  equilibrium 
(see  Fig.  5).  Methane  shows  a  strong  temperature  influence  at 
normal  operating  conditions.  The  concentration  doubles  from 
10  to  20%  if  the  temperature  decreases  about  70  K,  i.e.  small 
deviations  of  the  reactor  temperature  cause  obvious  changes 
in  the  methane  concentration.  Precise  agreement  of  simu¬ 
lated  and  experimental  CH4 -concentration  therefore  requires 
a  precise  temperature  simulation.  For  online  application  this 
can  be  achieved  by  an  observer,  which  combines  the  model 
with  temperature  sensor  data. 

The  dynamic  behaviour  for  both  operating  points  is  simi¬ 
lar: 

•  Increasing  the  natural  gas  flow  rate  does  hardly  influence 
the  thermal  system  but  causes  an  increased  outlet  methane 
concentration. 

•  On  the  other  hand  increasing  the  water  flow  rate  leads 
to  lower  evaporator  and  also  reactor  temperatures  due  to 
the  higher  evaporation  heat.  Initially  the  outlet  hydrogen 
concentration  is  increased  and  the  methane  concentration 
decreased,  but  this  is  later  on  compensated  by  the  lower 
temperature. 

•  Finally  increasing  the  burner  temperature  results  in  a 
higher  reactor  temperature  and  therefore  increased  hydro¬ 
gen  and  decreased  methane  concentration. 


6.  Conclusion 

We  described  the  derivation  of  a  dynamic  reformer  model 
based  on  a  theoretical,  lumped  element  approach.  The  model 
structure  and  the  parameter  determination  were  explained. 
Experimental  results  of  a  demonstration  plant  were  shown 
and  compared  with  simulation  results. 

The  presented  data  shows  the  good  performance  of  our 
model  despite  its  simple  structure.  Temperatures  and  gas 
compositions  agree  over  a  wide  operation  area  including  plant 
startup.  It  is  also  shown  that  for  the  outlet  gas  composition 
the  assumption  of  chemical  equilibrium  holds. 
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